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The Gamow- Teller response is astrophysically important for a number of nuclides, particularly 
around iron. The random phase approximation (RPA) is an efficient way to generate strength 
(N. distributions. In order to better understand both theoretical systematics and uncertainties, we 

{^jqj compare the Gamow- Teller strength distributions for a suite of nuclides and for a suite of interactions, 

r-{ _ including semi-realistic interactions in the lp-Of space with the RPA and a separable multi-shell 

. interaction in the quasi-particle RPA. We also compare with experimental results for GT_ on 54Fe. 

'. I. INTRODUCTION 

■ Gamow- Teller (GT) electron-capture transitions, caused by the ar± operator, are some of the most important 
nuclear weak processes in astrophysics. For a review of spin-isospin transitions see Ref. [l|. The GT transitions in 
/p-shell nuclei play important roles at the core collapse stages of supernovae, specially in neutrino induced processes. 
One of the factors controlling the gravitational core-collapse of massive stars is the lepton fraction; the lepton fraction 
in turn is governed by beta-decay and electron-capture rates among iron-regime nuclides. A primary and non-trivial 
contribution to the weak rates is the distribution of GT strength. GT strengths have important implications in other 

. astrophysical scenarios as well, such as explosive nucleosynthesis in O-Ne-Mg white dwarfs (see Ref. [2| and references 

■ therein) . 

While GT distributions have been extracted experimentally (but not in a model- independent manner), astrophysical 
CO • calculations rely either upon crude estimates or upon more detailed microscopic calculations. The main difficulty 
with both experiment and theory is that the strength distribution connects to many states. Further in astrophysical 
environments one needs finite temperature GT strength functions as the temperature is high enough for excited states 
OO 1 in the parent to be thermally populated. 

The isovector response of nuclei may be studied using the nucleon charge-exchange reactions {p, n) or (n,p); by 
other reactions such as ( 3 He,i), (<i, 2 He) or through heavy ion reactions. The 0° GT cross sections (AT = 1, AS = 
1 , 1,AL = 0, Ohuj excitations) are proportional to the analogous beta-decay strengths. Charge-exchange reactions at 
small momentum transfer can therefore be used to study beta-decay strength distributions when beta-decay is not 
energetically possible. The (p, ri), ( 3 He,t) reactions probe the GT_ strength (corresponding to /3~-decay) and the 
(n,p), (d, 2 He) reactions give the strength for /3 + -decay/electron capture, i.e. GT + strength. The study of (p, n) 
reactions has the advantage over /3-decay measurements in that the GT_ strength can be investigated over a large 
region of excitation energy in the residual nucleus. On the other hand the (n,p) reactions populates only T = T + 1 
states in all nuclei heavier than 3 He. This means that other final states (including the isobaric analog resonance) are 
forbidden and GT + transitions can be observed relatively free of background. The study of these reactions suggest 
that a reduction in the amount of GT strength is observed relative to theoretical calculations. The GT quenching is 
on the order of 30-40 % 0. 

Theory for GT transitions falls generally into three camps: simple independent-particle models (e.g. Ref. [![); 
full-scale interacting shell-model calculations; and, in between, the random-phase approximation (RPA) and quasi- 
particle random-phase approximation (QRPA). Independent-particle models underestimate the the total GT strength, 
because the Fermi surface is insufficiently fragmented, while also placing the centroid of the GT strength too high for 
even-even parent nuclides and too low on odd-A and odd-odd parents [5|. Full interacting shell- model calculations 
are computationally demanding, although one can exploit the Lanczos algorithm, commonly used in large shell-model 
diagonalization Q , to efficiently generate the strength distribution Q ; for medium-mass nuclei one still needs to choose 
from among a number of competing semi- realistic/semi-empirical interactions. RPA and QRPA can be thought of as 
approximations to a full shell-model calculation and are much less demanding computationally. 

In this paper we compare GT strength distributions for a suite of iron-region nuclides relevant to astrophysics: 
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54 ' 55,56 Fe, and 56 > 58 Ni. For each of these nuclides we compute the GT strengths in several RPA scenarios. Each RPA 
scenario is in occupation rather than configuration space, which is appropriate inasmuch as the GT operator only 
affects spin and isospin. (In fact each scenario properly speaking is proton- neutron RPA or QRPA, as the RPA/QRPA 
phonon operators change protons into neutrons or vice- versa.) The scenarios, which will be described in greater detail 
in subsequent sections, are : 

• pn-RPA in a major harmonic oscillator shell, that is, the lp-Of shell, with three different semi-realistic/semi-empirical 
interactions [sj-TToj . For details see Section II. 

• pn-QRPA in a multi-shell single-particle space with a schematic interaction that has been previously applied to 
similar calculations [ll|. For further details we refer to Section III. 

These particular scenarios were chosen because of the availability of codes; one could imagine a larger set of scenarios 
(e.g., pn-QRPA with semi- realistic interactions) but relevant codes either do not exist or are not available to us. 

These calculations will help to understand systematic similarities and differences between (a) different lp-Of shell- 
model interactions and (b) between OTiU shell-model calculations against multi-shell calculations with a separable 
interaction. For example, for some cases in the ls-Od shell the separable interactions yield a larger total strength and 
a higher centroid [13] than shell-model calculations. 

Section IV presents the results and discussions on use of various pn-RPA schemes. We finally present the summary 
and conclusions in Section V. 



II. THE RANDOM PHASE APPROXIMATION WITH SHELL-MODEL INTERACTIONS 

The configuration-interaction (CI) shell model solves the many-body problem in a large basis of Slater determinants 
using the occupation representation. One advantage of the CI shell model is that it can use arbitrary two-body (or 
even higher-order) interactions and gives explicit wavefunctions for excited states as well as the ground state. In 
addition, because the GT operator is err and does not affect coordinate space wavefunctions, the CI shell model is 
well suited for GT transitions. The drawback of the CI shell model is that even with including just a few, or even one, 
harmonic oscillator shell, the basis dimensions can be huge (10 9 or greater), making such calculations computationally 
intensive. 

The Hamiltonian for the shell-model is written in occupation space 0, [H, HH , i-e., 

H = ^ e a n a + i ^2 V abc dc\c\cdCd (I) 

a abed 

where the creation and annihilation operators c a ,c a represent single-particle states with good angular momentum, 
and where e Q are the single particle energies and V a bcd are the two-body matrix elements. 

It is possible to solve the RPA matrix equation in a shell-model representation, using the single-particle energies 
and two-body matrix elements above. A recent series of papers showed that RPA is a reasonable, if not perfect, 
approximation to the numerically exact results, comparing ground state correlation energies fl5j ] and charge-conserving 
fl6| and charge-changing [I?) transitions; in the last case it was found that allowing the Hartree-Fock state to be 
deformed improved pn-RPA calculations of GT strength distributions. 

The first step is a Hartree-Fock calculation, which introduces a unitary transformation on the single-particle states, 



4 = I>««4 (2) 



These states are divided into occupied (hole) states, labeled by m, and unoccupied (particle) states labeled by i, and 
the transformation matrix D is chosen such that the energy of the Slater determinant |HF) = rim^mlO) minimizes 
the energy (HF|ff |HF). 

With the Hartree-Fock solution in hand, one finds excited states (and the correlation energy in the ground state, 
although that does not concern us here) by treatin g th e energy surface in the vicinity of the Hartree-Fock state as 
quadratic. This leads to the RPA matrix equations |13j . For charge-changing interactions such as Gamow- Teller, the 
RPA matrix equations take the form: 



j^np,pn Q 
j\pn,np 

_fipn,np q 




/ X(pn) \ / X{pn) \ 

X(np) X(np) 

Y{np) Y(np) ' [6) 

\ Y(pn) J \ Y(pn) J 
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where the definitions for A pn ' np and B np,pn matrices are similar to the regular proton-neutron conserving formalism, 
where one approximates |RPA) ps |HF): 

^™ = <HF|[^ m , [if,4^]]|HF> = {el - e?)5 mn 5 l3 - VZ, t , (4) 

B%% = -(HFlf^TT,, \^y h U\\\m) = -Vf n n jm . (5) 

The matrices A pn,np and B pnjip are defined similarly, but are distinct unless Z = N; in fact, they have different 
dimensions unless Z — N. Let N£, N% be number of proton particle and hole states, respectively, and Ng, N£ the 
number of neutron particle and hole states. Thus the vectors X(pn) and Y(np) are of length N^N^ while vectors 
X(np), Y(pn) are of length N"N^\ the two lengths are unequal unless Z — N . Similarly, A np,pn is a square matrix 
of dimension N^N^ while A pn ' np is a square matrix of dimension NpNfi, while B np,p ™ is a rectangular matrix of 

dimension N%N£ x N»N%, and B pn - np = (B np ' pn f . 
The transition strength is given by 

(RPAIOlA^ii.jv^i)) = (RPA|[0,/?t]|RPA> = ]T (X^(pn/np)O m i + Y^{pn/np)O im ) . (6) 

mi 

For more details consult [13] • 

For this paper we use three different semi- realistic/semi-empirical shell model interactions. All three interactions 
started from realistic nucleon-nucleon interactions, from which an effective interaction (e.g., a G- matrix) was derived. 
At this point the interaction is expressed numerically as two-body matrix elements Vjxia-b, cd). The interactions 
were all then further modified in order to fit experimental spectra; as is well-known to the shell-model community, 
most of the modification were to the "monopole" parts of the interaction, which are related to properties of the 
mean-field. All three interactions are similar, but have different starting points and were fitted to different data sets, 
with the following semi- realistic/semi-empirical interactions: the modified Kuo-Brown interaction KB3G [8| and the 
Brown-Richter interaction interaction FPD6G Q and the Tokyo interaction GXPF1 [l(|; the names do not signify 
much except that PF/FP refer to the fp shell. 

III. THE QUASI-PARTICLE RANDOM PHASE APPROXIMATION WITH A SEPARABLE 

INTERACTION 

For an alternate approach, we used the quasi-particle proton-neutron random phase approximation (pn-QRPA) 
with a separable interaction of the form 

fjQRPA = jj S p + ypa lr + yph + ypp^ (7) 

where H sp is the single-particle Hamiltonian, V pair is the pairing force, Vq T and are the particle-hole (ph) and 
particle-particle (pp) components, respectively, of the GT force (<?f) . We diagonalized our Hamiltonian in three 
consecutive steps as outlined below. 

Single-particle energies and wave functions were calculated in the Nilsson model which takes into account nuclear 
deformation [l8| . The transformation from the spherical basis to the axial-symmetric deformed basis can be written 
as [HI 

dl a = Yl D T ad U' ( 8 ) 

3 

where eft and $ are particle creation operators in the deformed and spherical basis, respectively; the transformation 
matrices D™ a were determined by diagonalization of the Nilsson Hamiltonian, and a represents additional quantum 
numbers, except m, which specify the Nilsson eigenstates. 

Pairing was treated in the BCS approximation, where a constant pairing force with the force strength G (G p and 
G n for protons and neutrons, respectively) was applied, 

ypair = _ G £ (-l) ,+ ^" 1 cj m C ] t 1 „(-l) , '+^-»'c / _ m ,C,„ l -, (9) 
jmj'm' 
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where the sum over m and m! was restricted to m,m' > 0, and I represents the orbital angular momentum. 
The BCS calculation gave the quasiparticle energies e ma . A quasiparticle basis was introduced via 



where m is the time-reversed state of m, and a^(a) are the quasiparticle creation (annihilation) operators which 
enter the RPA equation. The occupation amplitudes u and v satisfy the condition u 2 + v = 1 and were determined 
by the BCS equations (see for example [l3j . page 230). 

In the pn-QRPA, charge-changing transitions are expressed in terms of phonon creation, with the QRPA phonons 
defined by 

bHn) = £(xn/i)fit4 - YS n {n)a n a ¥ ). (12) 

pn 

The sum in Eq. (j!2[) runs over all proton- neutron pairs with /x = m p — m n = -1, 0, 1, where m p ( n ) denotes the third 
component of the angular momentum. The ground state of the theory is dchncd as the vacuum with respect to the 
QRPA phonons, bui(fi)\QRPA >= 0. The forward- and backward-going amplitudes X and Y are eigenfunctions of 
the RPA matrix equation 



A 


B ' 




' X ' 




' X ' 


-B 


A 




Y 




Y 



(13) 



where to are energy eigenvalues of the eigenstates and elements of the two submatrices are given by 

= 5(pn,pn)(c p + e n ) + V p P ny n '( u p u n u p /U n' + v p v n v p 'V n/ ) + 

), (14) 



B pn ,p>n> = Vpn, p >n'( U P U n v p' v n' + V p V n U p ,U n ') - V^ p , n , (u p V n V p > U n < + V p U n U p 'V n >). (15) 

The backward-going amplitude Y accounts for the ground-state correlations. It is essential to note however, that 
the derivation of the QRPA matrix requires ground-state correlations to be only a small correction. It should be noted 
that \Y\ ^ |A| does not imply that ground-state correlations are negligible, since for the calculation of /? transition 
matrix elements one always must consider products of the form uvY and u'v'X. Especially in f3 + decay, uv can be 
larger than u'v'\ thus ground-state correlations cannot be neglected. The RPA equation is constructed and solved 
for each value of the projection /i, i.e., /j, = -1, and +1. The equation gives identical eigenvalue spectra for /i = -1 
and /i = +1, and eigenvalues for /i = are always two-fold degenerate, because of the axial symmetry of the Nilsson 
potential. (Hereafter, /i will be suppressed if not otherwise stated, since the following formulae hold for each /j,). 

In the pn-QRPA formalism proton-neutron residual interactions occur in two different forms, namely as particle- 
hole (ph) and particle-particle (pp) interaction. Both the particle-hole and particle-particle interaction can be given 
a separable form. 

In the present work, in addition to the well known particle-hole force pjl, [2(| 

t^ = 2x£(-ir^, (16) 

with 

% = 0p m pl*-°'A'bn m n) 5 im p 6 inm„, (17) 

j p m p j„m„ 

the particle-particle interaction, approximated by the separable force [U [22j 

t^ = -2 K £(-l)^P^, (18) 
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with 

Pt= E (inm n |(f-^)+|i P m p )(-l) i »+^- m "ct pmp 4 n _ mn) (19) 

j P m p j n m n 

was taken into account. The interaction constants x an d K in units of MeV were both taken to be positive. The 
different signs of V pp and V ph reflect a well-known feature of the nucleon-nucleon interaction; namely, that the ph 
force is repulsive while the pp force is attractive. The selections of these two constants were done in an optimal 
fashion. For details of the fine tuning of the GT strength parameters, we refer to [23|, |24j. The values of x an d K, 
along with the value of deformation parameter, used in the current work, are shown in Table HI 

Using a separable interaction implies that the pn-QRPA calculations can be solved in a much larger single-particle 
basis than with a general/semi-realistic interaction; in this case we used up to 7 fiuj shells. Such calculations have 
been used extensively in computing GT transitions for astrophysical applications for a wide variety of nuclide (e.g. 

Matrix elements of the forces which appear in RPA equation (|14|) .([T5 |) are separable, 

<W = +2x/p»M/ P '»'(M), (20) 

Klp'n' = -2«/j»M/y B 'W> (21) 

with 

fpniv) = E D 7p* a * '-D™""" (i P TOp|t-o- M |j n m n ), (22) 

which are single-particle GT transition amplitudes defined in the Nilsson basis. For the case of separable forces, the 
matrix equation (|13[) reduces to an algebraic equation of second order (when k = 0) and with a finite value of k it 
transforms to a fourth order equation. Methods of finding roots of these equations can be seen in Ref. [1JJ. F° r 
details on QRPA model parameters we refer to [28j . 

The purpose of this paper is to compare general trends of these calculations with the Ohuj calculations using a more 
general, realistic interaction described in the previous section. 



IV. RESULTS AND COMPARISON 



Using the Gamow- Teller operator or± yields the Ikeda sum rule [29( for a parent nucleus with Z protons and N 
neutrons: 

E#(GT_) -^B{GT+) = 3(N - Z), (23) 

in our calculations. For use in astrophysical reaction rates, and to compare to experimental data, we have to multiply 
our calculated strengths by (gA/gv) 2 = (— 1.254) 2 [3(|, and then by an additional quenching factor of 0.6 Q typical 
for nuclei. 

The ultimate goal is to provide reliable weak rates for astrophysical environments, many of which cannot be 
measured experimentally. Even theoretically this is a complex and difficult issue. For example, /3-decay and capture 
rates are vexponentially sensitive to the location of GT + resonance while the total GT strength affect the stellar rates 



in a more or less linear fashion 31]. In sufficiently hot astrophysical environments one must include rates with an 
excited parent state. But rates off excited states are difficult to get: an (n,p) experiment on a nucleus (Z,A) shows 
where in [Z — 1, A) the GT + centroid corresponding only to the ground state of (Z,A) resides. The calculations 
described in this paper are also limited only to ground state parents, although we hope to tackle excited parents in 
the future. 

For this paper we focus on the variation in Gamow- Teller strengths from different RPA calculations, to give us an 
idea of the theoretical uncertainty. 

Table |TI] shows the mutual comparison of the various RPA models used in this project. This table shows the values 
of the centroids, widths and total strength values of the calculated GT distributions, both in j3~ and electron capture 
directions, for various iron-regime nuclei. It is clear from Table |TT] that the GXPF1 interaction calculates the biggest 
value of the total GT strength (unquenched). The lowest total strength is calculated by the pn-QRPA model except 
for the case of 58 Ni where the FPD6 interaction bags the lowest total. Regarding the calculation of centroids in various 
RPA models, we note that the pn-QRPA calculated centroid resides at lower energy in daughter, except for the case 
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of 55 Fe where the KB3G interaction calculates the lowest centroids in both directions. On the other hand the GXPF1 
calculates the highest centroids except for the case of 55 Fe where the pn-QRPA model tops the chart. One also notes 
that the pn-QRPA calculated GT strength distributions tend to have a larger width, specially for the case of 55 Fe in 
the j3~ direction (of around 15 MeV). The calculated GT strength distributions using different interactions will next 
be discussed below. 

An obvious question would be how the various RPA scenarios compare with the measured data as well as full shell- 
model diagonalization. Thus we also show in TableUTlhow the values of calculated GT + centroid and total GT strengths 
(both in P~ and electron capture directions), using various RPA models, compare with the available experimental 
data. For the sake of comparison we also include the shell model calculation using the GXPF1J Hamiltonian [32[ taken 
from Table I of Ref. [33[ (the centroid value was not available). The authors claimed that the GXPF1J interaction 
leads to spreading of calculated strength and better reproduction of observed strength in Fe and Ni isotopes. For the 
f3~ side the measured data for 54 Fe, 56 Fe and 58 Ni were taken from Ref. [34|. We note that a recent high-resolution 
( 3 He, t) charge- exchange reaction on 54 Fe performed by Adachi and collaborators [35[ report a much lower value of 
Y A B(GT-) = 4.00 ± 0.37 up to 12 MeV in 54 Co. For the (3 + side, experimental data for 54 Fe were taken from Ref. 
[3a ] while that for 56 Fe and 58 Ni were taken from Ref. f37|. It is to be noted that we used a quenching factor of 0.6 
for the calculated GT strength using the pn-QRPA model [38[ which is normally done in stellar weak rate calculations 
and also discussed in Section I. Note that the shell model interactions [H[ calculated strengths were quenched by a 
universal quenching factor of 0.55 rather than 0.6. Table UT1 shows that the pn-QRPA model calculates the centroid at 
a much lower energy than other shell model interactions. The comparison is exceptionally good for the case of 58 Ni. 
The GXPF1 interaction calculates the highest centroid in daughter nuclei. The shell model interactions calculate 
much better total strengths in comparison with measured values. One should also keep in mind the uncertainties 
present in measurements where various energy cutoffs are used as a reasonable upper limit on the energy at which GT 
strength could be reliably related to measured AL = cross-sections as well as slightly different values of quenching 
factor used in pn-QRPA and shell model interactions before comparing the calculated numbers with experimental 
data. The model independent Ikeda sum rule is satisfied in all models (except for the odd-A case in the pn-QRPA 
model). 

Fig. [U and Fig. [5] show the calculated GT strength distribution of 54 Fe in various RPA scenarios. Fig. [1] displays 
the calculated GT strength distribution in the electron capture direction whereas Fig. [2] shows similar calculation in 
the j3~ direction. In the inset of Fig. [5] we also show the recently measured high-resolution ( 3 He, t) data on 54 Fe by 
Adachi and collaborators [35| . Fig. [5] sh ows the cumulative strength distributions, in the /3 _ direction, for various 
RPA calculations and measured data [35(. Fig. [3] depicts the mutual comparison of various calculations with measured 
data in a better fashion. It is noted that the pn-QRPA model calculates GT transitions at low excitation energy. It 
is further noted that the QRPA calculated distribution is better fragmented than other RPA models and follows the 
trend of the measured data, albeit with a much higher magnitude of strength distribution. 

For the remaining cases we decided to show only the cumulative GT strength distributions to save space. Fig. 0] 
display the cumulative strength distributions in a two-panel frame. The upper panel shows the cumulative strength 
distribution in the f3~ direction using the different RPA models whereas the lower panel displays the results in the 
electron capture direction for the case of 55 Fe. The pn-QRPA calculates high-lying transitions in 55 Co (upper panel). 
Whereas no significant strength is calculated between 10 - 35 MeV, around 6 units of strength is concentrated in the 
energy region 35 - 44 MeV in daughter nucleus. Accordingly the centroid of this distribution, as calculated by the 
pn-QRPA model, is around 23 MeV in 55 Co (see Table ITT|) . The KB3G interaction saturates first to its maximum 
strength in both directions and gives the lowest values for the corresponding centroids of the GT distribution function. 

A similar comparison of cumulative strength functions for the case of 56 Fe is shown in Fig. [5j The figure reveals 
that the pn-QRPA model peaks at a faster pace compared to other RPA models. Correspondingly the pn-QRPA 
calculates the lowest values for the centroids in both directions for 56 Fe (see Table HI)) . It can be seen that the GT 
strength resides at much higher energies in the daughter nuclei in the j3~ direction for all models. 

Fig. [6] shows the GT strength distribution of N = Z nucleus 56 Ni. Here one notes that for the models used, 
all strength resides at one energy level in daughter nucleus for the FPD6, GXPF1 and KB3G interactions. This 
energy is 10.6, 11.5 and 9.8 MeV for the FPD6, GXPF1 and KB3G interactions, respectively. The recent shell model 
calculation by Suzuki et al. [Hj], using the GXPF1J interaction, showed that the GT strength distribution in 56 Ni 
is more fragmented and different from previous shell model calculations resulting in enhanced production yields of 
heavy elements. The authors further commented that the calculated GT strength using the GXPF1J interaction was 
found to be more fragmented with a remaining tail in the high excitation energy compared with that obtained by the 
KB3G interaction. For the case of 56 Ni their calculated strength was fragmented into two peaks whereas the total 
calculated strength was 11.32 (unquenched) . On the other hand we note from Fig. [6] that the pn-QRPA calculated 
strength is still more fragmented over a range of energies with two distinct peaks at 5.7 MeV and 10.6 MeV in the 
daughter nucleus. The consequences of calculating a much more fragmented strength distribution for 56 Ni, using the 
pn-QRPA model, may also have interesting scenario for heavy element nucleosynthesis and requires further attention 
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in nuclear network calculations. 

The cumulative strength distributions for the case of 58 Ni is shown finally in Fig. [7J Here one can see that the 
FPD6, GXPF1 and KB3G models calculate strength up to higher energies in 58 Cu as compared to the pn-QRPA 
model (upper panel). Consequently the pn-QRPA models calculate the GT_ centroid at a much lower energy of 
around 5 MeV in 58 Cu. The GXPF1 calculated centroid is thrice this value. A similar situation is seen for the GT 
strength distribution in the electron capture direction (lower panel). Once again the pn-QRPA model calculates the 
lowest whereas the GXPF1 model calculates the highest values for the GT + centroid. It is the case of 58 Ni where one 
sees the largest difference in the placement of centroid using the QRPA and various RPA interactions. 

V. SUMMARY AND CONCLUSIONS 

Under astrophysical conditions, both the electron capture and beta decay of fp-shell nuclei depend heavily on the 
centroid placement and total strength of the calculated Gamow- Teller strength distributions. In this work we presented 
a comparative study of the Gamow- Teller strength distributions for a suite of astrophyiscally important fp-shell nuclide 
( 54 ' 55 ' 56 Fe, and 56 > 58 Ni) using a suite of interactions, including semi-realistic interactions in the lp-Of space with the 
RPA and a separable multi-shell interaction in the quasi-particle RPA. Where possible, we also presented comparison 
with measured data. We further compared and contrasted the statistics of calculated GT strength functions using 
various pn-RPA schemes in this paper. Our calculations satisfied the model independent Ikeda sum rule. Work is 
currently in progress for other important odd-A and odd-odd cases. 

The QRPA model places the centroid at much lower energies in daughter nuclei as compared to other RPA inter- 
actions. This tendency of QRPA model favors higher values of electron capture rates in stellar environment. On the 
other extreme, the GXPF1 interaction usually leads to placement of GT centroid at much higher energies in daughter 
compared to other pn-RPA interactions. 

The present study showed that the total strengths, using various RPA interactions, were in better agreement with 
the measured data when compared to the QRPA calculated strength. Further the width of the strength functions 
calculated within the QRPA scenario was much larger than those calculated with other RPA interactions. For the 
special N — Z nucleus 56 Ni, the QRPA model calculated Gamow- Teller strength function was well fragmented as 
compared to other RPA interactions (including the recently used GXPF1J interaction) and may lead to interesting 
consequences for heavy element nucleosynthesis. 
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TABLE I: Values of GT strength force constants and deformation parameters used in the pn-QRPA calculation 



Nucleus 


X 


K 


S 


54 Fe 


0.6500 


0.0775 


0.195 


55 Fe 


0.6500 


0.0775 


0.083 


56 Fe 


0.6500 


0.0775 


0.239 


56 Ni 


0.5000 


0.0650 


0.011 


58 Ni 


0.5000 


0.0650 


0.183 



TABLE II: Statistics of calculated GT strength distributions for nuclei using different RPA models given in second column. 
The GT strength values for shell-model calculations using GXPF1J interaction were taken from Table I of Ref. [33J, while 
those for experiment are from [341 ]. 



Nucleus 


Model 


E(GT + ) 


Width(GT + ) 


V B(GT + ) 


E(GT_ 


) Width(GT_ 


) Vfl(CT_) 






MeV 


McV 


arb. units 


McV 


MeV 


arb. units 


54 Fe 


















FPD6 


7.50 


1.01 


6.66 


11.92 


2.56 


12.32 




GXPF1 


8.42 


1.00 


8.16 


12.97 


2.59 


13.82 




KB3G 


7.13 


0.83 


7.27 


11.02 


2.34 


12.95 




QRPA 


5.80 


3.38 


5.65 


8.03 


4.40 


11.65 




GXPF1J (SM) 






4.0 






7.3 




EXP 


3.7 ± 0.2 




3.5 ± 0.7 






7.8 ± 1.9 


55 Fe 


















FPD6 


4.42 


1.43 


5.05 


14.87 


3.02 


13.54 




GXPF1 


5.07 


1.26 


7.11 


15.93 


3.24 


15.60 




KB3G 


4.03 


1.08 


6.01 


13.96 


2.95 


14.50 




QRPA 


7.12 


1.75 


4.69 


23.07 


15.23 


13.59 


56 Fe 


















FPD6 


6.12 


1.75 


3.74 


12.79 


3.00 


15.07 




GXPF1 


6.73 


1.38 


5.82 


14.09 


3.39 


17.14 




KB3G 


5.49 


1.32 


4.88 


12.54 


3.13 


16.20 




QRPA 


3.14 


1.53 


3.71 


7.79 


3.79 


15.71 




GXPF1J (SM) 






2.9 






9.5 




EXP 


2.6 ± 0.2 




2.9 ± 0.3 






9.9 ± 2.4 


56 Ni 


















FPD6 


10.62 


0.00 


9.97 


10.62 


0.00 


9.97 




GXPF1 


11.54 


0.00 


11.92 


11.54 


0.00 


11.92 




KB3G 


9.77 


0.00 


10.71 


9.77 


0.00 


10.71 




QRPA 


6.32 


1.91 


8.87 


6.32 


1.91 


8.87 


58 Ni 


















FPD6 


6.15 


0.92 


6.86 


13.94 


2.49 


12.52 




GXPF1 


6.79 


0.58 


9.58 


14.99 


3.34 


15.25 




KB3G 


5.02 


0.70 


7.72 


13.97 


2.55 


13.38 




QRPA 


3.57 


1.91 


7.82 


4.97 


2.82 


13.82 




GXPFIJ(SM) 






4.7 






8.0 




EXP 


3.6 ± 0.2 




3.8 ± 0.4 






7.4 ± 1.8 
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FIG. 1: Calculated GT strength distributions in the electron capture direction for Fe using different pn-RPA scenarios. For 
details of interactions used see text. The abscissa represents energy in daughter. 
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FIG. 2: Same as in Fig. [T]but for calculated GT strength distributions in the /3 -decay direction. In inset the recent measured 
data of Ref. [35[ is shown. 
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FIG. 3: Cumulative GT strength distributions for 5 Fe using different pn-RPA scenarios in the ft -decay direction. For details 
of interactions used see text. Experimental data (exp) is taken from Ref. [33|. The abscissa represents energy in daughter. 
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FIG. 4: Cumulative GT strength distributions for 55 Fe using different pn-RPA scenarios. For details of interactions used see 
text. The abscissa represents energy in daughter. 
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FIG. 5: Same as in Fig. [J but for 56 Fe 
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FIG. 6: Same as in Pig. H but for the N = Z nucleus 56 Ni 
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FIG. 7: Same as in Fig. g] but for 58 Ni 




